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FAR  FIELD  COMPUTATIONAL  BOUNDARY  CONDITIONS 


FOR  INTERNAL  FLOW  PROBLEMS 
A.  VERHOFF 

MCDONNELL  AIRCRAFT  COMPANY 
ST.  LOUIS,  MISSOURI 

ABSTRACT 

Far  field  computati onal  boundary  conditions  for  2D  internal  flow  problems 
are  developed  from  analytic  solutions  of  the  linearized  Euler  equations.  The 
Euler  equations  are  linearized  about  a  constant  pressure,  rectilinear  f 1 ow 
which  nay  have  streamwi se-normal  variations  in  temperature  and  velocity  as  a 
result  of  entropy  production  in  the  nonlinear  computational  region.  The  bound¬ 
ary  procedure  can  be  used  with  any  numerical  Euler  solution  method  and  allows 
computational  boundaries  to  be  placed  much  closer  to  the  nonlinear 
region  of  interest. 


1.  INTRODUCTION 


Numerical  solution  procedures  for  nonlinear  fluid  dynamic  equations 
usually  use  one  or  more  artificial  computational  boundaries  located  at  some 
distance  from  the  primary  region  of  interest  in  order  to  limit  the  physical 
domain  to  finite  size.  If  the  flow  crossing  such  a  boundary  (either  inflow  or 
outflow)  is  subsonic,  then  some  type  of  computational  boundary  conditions  must 
be  imposed  which  simulate  the  influence  of  the  true  far  field  conditions  at 
infinity.  These  boundary  conditions  must  be  such  that  waves  crossing  the 
boundary  do  not  produce  erroneous  reflections  back  into  the  computational 
field  to  degrade  the  calculations.  It  is  generally  acknowledged  that  simply 
imposing  free  stream  conditions  (or  conditions  at  infinity)  at  computational 
boundaries  is  usually  inappropriate.  Standard  practice  has  consisted  of 
locating  the  boundaries  quite  far  from  the  region  of  interest  in  an  attempt  to 
simplify  the  boundary  condition  models  and  minimize  any  effects  of 
inconsistent  modeling.  The  net  effect  is  a  significant  increase  in  the  number 
of  grid  points  required  for  an  accurate  flowfield  calculation. 

A  boundary  modeling  procedure  for  two-dimensional  internal  flows  is 
presented  which  alleviates  the  difficulties  mentioned  above  and  also  allows 
the  computational  boundary  to  be  located  much  closer  to  the  nonlinear  region 
of  interest.  The  procedure  is  limited  to  steady,  inviscid  flow,  although  the 
flow  can  be  rotational.  It  represents  a  logical  extension  of  the  so-called 
character i sti c  boundary  conditions  commonly  used  with  inviscid  numerical 
solution  methods.  Extension  to  axisymmetric  or  three-dimensional  flows  is 
straightforward. 

The  analysis  presented  here  is  based  on  the  Riemann  variable  formulation 
of  the  Euler  equations  given  in  Reference  1.  This  represents  a  natural 
starting  point  because  the  characteri sti c  (or  zero-order)  boundary  conditions 
mentioned  above  are  expressed  in  terms  of  Riemann  variables.  The  equations 
are  linearized  about  a  constant  pressure,  rectilinear  flow  condition,  which 
truly  represents  conditions  at  infinity.  These  linearized  equations  are 
assumed  applicable  in  the  far  field  region  beyond  a  computational  boundary. 
Within  the  nonlinear  computational  domain,  strong  entropy-producing  (i.e., 
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rotational)  effects  can  exist  which  create  variations  in  density,  velocity, 
etc.  in  the  far  field  in  the  streamline-normal  direction  which  are  not 
necessarily  small  perturbations.  Such  variations  are  modeled  in  the  present 
analysis . 

The  linearized  equations  are  solved  analytically  using  Fourier  analysis 
techniques  as  outlined  in  Reference  2.  These  solutions  are  coupled  to  the 
nonlinear  numerical  solution  to  provide  a  smooth  transition  across  the 
boundary  to  the  true  far  field  conditions  at  infinity.  The  coupling  is 
accomplished  by  the  boundary  conditions.  The  underlying  principle  is  that  the 
streamwise  variations  of  both  upstream  and  downstream  running  waves  should 
decay  to  zero  at  infinity.  These  first-order  boundary  conditions  provide 
distributions  of  flow  quantities  to  be  imposed  along  the  boundary,  not 
constant  conditions.  They  represent  a  logical  extension  of  the  zero-order  (or 
characteristic)  boundary  conditions.  Furthermore,  the  boundary  analysis  can 
be  coupled  with  any  inviscid  numerical  solution  method. 

The  boundary  condition  analysis  has  been  applied  to  two-dimensional 
duct  flow  and  to  cascade  flow,  where  conditions  are  periodic  in  the  direction 
of  the  blade  row.  For  duct  flow,  both  isentropic  and  non-isentropi c  boundary 
conditions  are  derived.  Only  isentropic  results  are  given  for  cascades. 
Extension  to  non-i sentropi c  cascade  flow  can  be  carried  out  by  following  the 
procedure  used  for  duct  flow. 

Numerical  results  are  presented  for  both  isentropic  and  non-isentropic 
duct  flow.  Results  obtained  using  the  first-order  boundary  conditions  are 
compared  with  those  using  the  zero-order  boundary  conditions.  Numerical 
results  for  cascade  flow  will  be  presented  in  Reference  3.  It  was  found  that 
the  size  of  the  computational  field  and  associated  number  of  grid  points 
needed  for  the  nonlinear  numerical  solution  could  be  reduced  significantly  by 
using  the  new  first-order  boundary  condition  procedure  with  no  loss  in 
numerical  solution  accuracy.  The  reduction  in  number  of  grid  points  was  as 
much  as  50  percent  in  some  cases.  The  additional  computational  effort 
required  by  the  new  boundary  procedure  is  small  (less  than  10  percent)  so  that 
a  significant  saving  in  overall  computational  effort  was  realized.  A  large 
portion  of  the  gain  is  due  to  the  sizeable  reduction  in  the  physical  extent  of 
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the  computational  field  and  the  fact  that  fewer  solution  iterations  are 
required  for  information  to  propagate  between  the  upstream  and  downstream 
computational  boundaries  leading  to  more  rapid  solution  convergence. 
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II  .  PERTURBATION  EQUATION  DERIVATION 


The  system  of  two-dimensional,  steady,  linearized  Euler  equations  which 
describe  first-order  spatial  perturbations  from  a  constant  pressure  state  will 
be  derived  in  this  section.  A  Riemann  variable  formulation  taken  from 
Reference  1  will  be  used  because  of  its  close  relationship  with  the  character¬ 
istic  (or  zero-order)  boundary  conditions  commonly  used  in  numerical  solution 
of  the  nonlinear  Euler  equations. 


The  two-dimensional  form  of  the  Euler  equations  is  (Reference  1) 


a  (S  - 


.  _L_  ill  .  zzl  aaS 
y-r  l3s  Y-l  3sJ  2  q  b  an 


aR  ,  .  aR  _  .  v-i  -  ,r 

H  (q  '  a)  al  2  a  (s 


ri9  +  . +  Ili.  Qa$ 

r-1  Las  Y-l  asJ  2  qab 


30 

3n 


30  ,  36 

at  q  as 


a£  aP 

yq  3n 


3S  .  n  3S 
at  q  as 


=  o 


(1) 

(2) 

(3) 

(4) 


Velocity  magnitude  and  speed  of  sound  are  denoted  by  q  and  a,  respectively, 
and  P  is  the  logarithm  of  pressure.  The  Extended  Riemann  Variables  Q  and  R 
are  defined  as 


Q  =  q  +  aS  (5) 

R  =  q  -  aS 

while  the  modified  entropy  is  defined  in  terms  of  pressure  p  and  density  p  as 

s  =  Trhr  [2y  " 109  (F>/p1f)]  (6) 

The  flow  angle  is  0,  time  is  denoted  by  t,  and  local  distances  along  and 
normal  to  the  streamline  direction  are  denoted  by  s  and  n,  respectively. 
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For  steady  flow  the  analysis  can  be  greatly  simplified  by  defining  a  new 
dependent  variable 

T  =  Q  -  R  (7) 

and  replacing  equations  (1)  and  (2)  by 

(M2  -  1)  +  (y-1)  Q  M  S  |£  =  0  (8) 

a2  +  ^  q2  =  1  (9) 

The  local  Mach  number  is  denoted  by  M.  Equation  (8)  is  obtained  by 
subtracting  equations  (1)  and  (2).  Equation  (9)  is  obtained  by  adding 
equations  (1)  and  (2)  and  integrating.  The  constant  of  integration,  which  is 
proportional  to  stagnation  temperature,  can  be  set  to  unity  by  proper  choice 
of  non-dimensional izing  quantities.  The  simplified  form  of  the  steady  Euler 
equations  is  then 


(MZ  '  1>  (y-1)  OM  Sf£=  0 


f^|2  98  +  — 2—  I  il  +  I  (c  -  2  ) 

n  9s  Y-1  T  an  S  v-W  An 


Y~1 y  9n 


ff=°  (12) 

a2  +  III  q2  =  !  (13) 

According  to  equation  (12)  entropy  remains  constant  along  streamlines. 

In  regions  of  the  flowfield  where  nonlinear  effects  are  weak,  the  flow  can  be 
treated  as  a  perturbation  to  a  constant  pressure,  rectilinear  flow.  Such 
regions  occur  near  and  beyond  far  field  computational  boundaries.  The  depend¬ 
ent  variables  in  equations  ( 10 ) - ( 12 )  can  then  be  expanded  in  asymptotic  series 

T  =  T  +  T.  +  T  +  ... 

“12 


S  =  $  +  S,  +  S,  +  ... 
“12 


e=e  +  e,  +  e,  +  . . 

“12 
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The  flow  direction  at  infinity  is  assumed  constant  and  denoted  by  0^;  the 
perturbation  quantities  T^,  and  8^  vanish  at  infinity.  Entropy  variation 
is  not  excluded  so  that  the  flow  can  be  rotational.  Furthermore,  entropy 
variations  can  be  strong  (i.e.,  not  small  perturbations)  so  that  and  are 
not  necessarily  constant,  but  may  vary  normal  to  the  streamline  direction. 

Note  that  T  depends  only  on  S  (because  p  is  constant). 

OO  00  00 

Consistent  with  (14),  spatial  derivatives  in  equations  ( 10 ) - ( 12 )  can  be 
approximated  by 

~  =  cos  ( l“t8j )■—  +  cos  0„  (t+6i)|y  +  ... 

(15) 

=  -cos  0M  ('t+0l)|^  +  COS  0oo  (1-T0l)|y  +  ... 

where  x  and  y  are  reference  Cartesian  coordinates,  e  is  measured  from  the  x 
axis,  and 

i  =  tan  0^  (16) 

If  expansions  (14)  and  (15)  are  introduced  into  equations  ( 10 ) - ( 12 ) ,  the 
resulting  first-order  perturbed  Euler  equations  are 


,M2  i \ rill  +  !Il 

^  J^3x  t  8y 


(t-1)  M«,  S„  [; 


^  «  OO  \  -t 

el  ST  '  ar)]  * 


301 


111  M2  C  rill  4.  iili  ,  1  f  aTl  ill  li  rilz 

2  Mco  [3x  T  av  ^  2a*  *  ay  T  ax  T,*,  *-ay 


T  3F])  + 


1  ,  r9Si  9Si  Si  as* 

y-l;  Lay  T  ax  J  s«  Lay 


-  '  aF^  =  0 


3Si  3Si  3Sa 

aF  +  T  W  +  81  [W 


as-  n 
T  aF]  =  0 
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Velocity,  speed  of  sound,  and  Mach  number  at  infinity  are  denoted  by  q  ,  a 

oo  oo 

and  Mm,  respectively,  and  may  vary  normal  to  the  streamline  direction.  The 
fact  that  all  dependent  quantities  are  convected  unchanged  along  rectilinear 
streamlines  at  infinity  has  been  incorporated  into  the  above  equations. 

Asymptotic  expansions  of  the  Riemann  variables  Q  and  R  can  also  be 
defined  as 


Q  =  Q.  +  Qi  +  Q2  +  • • • 

R  =  Rco  +  R1  +  R2  + 


(20) 


Using  the  definition  (7)  and  the  expansion  (14)  for  T,  it  follows  that 


(21) 


Introducing  the  expansions  (14)  and  (20)  into  the  algebraic  equation  (13) 
gives  the  first-order  relationship 

NL  (Qi-Ri)  +  ^  S«,  m£  (Qi+Rj)  -  2q«  Sj  =  0  (22) 

This  will  be  used  later  in  Section  IV  where  the  boundary  conditions  are 
derived. 
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Ill .  SOLUTION  OF  FIRST-ORDER  EQUATIONS 

Solutions  of  the  first-order  equations  ( 17 ) - ( 19 )  are  developed  in  this 
section  for  two-dimensional  duct  flows  confined  by  parallel  walls  and  for 
infinite  cascade  flows.  Duct  flow  solutions  are  obtained  for  both  isentropic 
and  non-isentropic  conditions.  The  cascade  flow  solution  is  provided  for 
isentropic  conditions. 

Duct  Flow 

For  two-dimensional  duct  flow  between  parallel  walls,  the  far  field  flow 
angle  6^  (and  consequently  t)  can  be  set  to  zero  by  aligning  the  coordinate  x 
with  the  duct  axis  (see  Figure  1).  The  width  of  the  duct  can  be  taken  as 
unity  without  loss  of  generality.  Equations  ( 17 ) - ( 19 )  then  reduce  to 

p  3Tl  dloo  901 

(M^  -1)  [~  +  ei  — ]  +  b-1)  Qm  :  =  0  (23) 


d  M  c  +  i  r^Il  li  .  a  !i  ^ 

2  q™  ”  “  9x  2  L9y  '  T^  9y  J  S„  9y 


+  111  a  (S  -  )  — —  =  o 

2  a“  y- 1  9y  U 


(24) 


3Si  aSoo 

— -  +  01  -  =  0 

9x  1  9y 


(25) 


Isentropic  Conditions  -  According  to  the  scaled  entropy  definition  (6),  the 

far  field  entropy  S  is  — r.  Equations  ( 23 ) - ( 25 )  then  simplify  for  isentropic 

<*>  y-1 

flow  to 


2qoX 


9x 


90  , 

W 

9Ti 


891  +  _1 _ 

9x  2q00Mw  9y 


=  0 


=  0 


(26) 

(27) 
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where 


(28) 


Note  that  q  and  M  are  constant. 

oc  oo 


Equations  (26)  and  (27)  can  be  solved  by  separation  of  variables  by 
assuming 


Tj  =  2  q^  eAx  F(y) 
0X  =  eXx  H(y) 


(29) 


where  X  is  the  unknown  separation  constant.  Equations  (26)  and  (27)  then 
reduce  to  the  ordinary  differential  equation  system 


p2  X  F  -  H1  =  0 
XH  +  F 1  =0 


(30) 


with  boundary  conditions 

H( 0 )  =  H ( 1 )  =  0  (31) 

Primes  denote  differentiation  with  respect  to  y.  Eliminating  F  gives  the 
single  equation 

h"  +  xV  H  =  0  (32) 

which  has  the  solution 

H  =  kj  sin  xpy  +  k ^  cos  xpy  (33) 

To  satisfy  the  boundary  conditions  (31),  the  constant  k^  must  be  zero  and 

X  =  t  *7  (n  =  1,  2,  ...)  (34) 
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Using  these  results,  the  general  solution  of  the  system  (26)  and  (27)  can 
be  written  as 


Tl 

el 


oo 


I  Kln 
1 


^gooMoo 

6 

si 


cos  nny 

n  nny 


rnix/p 

e  + 


Z  K2n 
1 


P 


cos  n-ny 


sin  nity 


-n-nx/ p 


(35) 


where  K^n  and  K^n  are  arbitrary  coefficients. 


Non-lsentropic  Conditions  -  Although  the  system  (23)-(25)  is  linear,  it  has 
non-constant  coefficients  since  is  independent  and  varies  with  the 
streamline-normal  direction  y.  The  thermodynamic  variables  are  related  by 

(36) 

where  is  the  constant  pressure  at  infinity. 

To  achieve  an  analytic  solution,  a  second  level  of  1 inearization  can  be 
introduced  by  defining  a  new  variable  o(y)  by 

s„  s  <l-°>  (3)) 

Since  a  is  generally  small,  equation  (36)  can  be  approximated  by 

111 

a«,  =  Poo  ^  (l+o)  (38) 
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An  approximation  for  q^  is  then  provided  by  equation  (9).  Reference 
quantities  defined  by 


a«  =  p0 


~  (l-o 

=  iJK 


can  be  introduced  to  simplify  the  notation.  The  approximate  thermodynamic 
relations  then  become 

a«,  =  a,,,  (1  +  o) 

q«  -  doo  (1  -  «  _  ) 

Y  1  M2 


rio  =  Me  (1  -  =f) 

Y  i  q*- 


Too  3oo 

Y-l 

All  quadratic  and  higher  terms  in  o  have  been  neglected. 


Introducing  the  above  approximations  (40)  into  the  governing  equations 
(23)-(25)  gives 


/„ 2  +  J4 _ £\  £Li  On  M  M  _ 2.\  -  n 

{p  Y-l  a2  3x  ’  2q"Mco  1  Y-l  q2J  ay  "  0 


On  M  /  1  4  "I  .  9T1  _  o'  9  ,  r  v 

2q«n„  (1  -  y_j  -2)  3x  +  3y  -  230,,  3y(oSj) 


Hi  .  _2_  do  _  Q 
8x  y-l  dy  * 
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The  parameter  g  is  defined  by  equation  (28)  with  M  substituted  for  M  . 

oo  oo 

Equation  (43)  suggests  that  is  of  higher  order  than  or  6^.  If  so,  then 
the  RHS  of  equation  (42)  can  be  set  to  zero  and  equation  (43)  becomes 
decoupled.  Validity  of  this  assumption  can  be  verified  later  by  evaluating 
numerical  results  obtained  from  this  analysis.  The  final  system  of  governing 
equations  is  therefore 


o .  ^1 


^  +  A  n>  >r  - ^  (i 


JL  -a,  =  o 

y-1  q2  ay 


-  -  4  n  301  3Ti 

2qJ4.  (1  -±1$  °  (45) 

OO 

After  these  equations  are  solved  for  T^  and  0^,  the  perturbation  variable 
can  be  obtained  from  equation  (43)  by  quadrature. 


As  for  the  isentropic  analysis  presented  above,  a  separation  of  variables 
solution  approach  is  again  appropriate.  Analagous  to  equations  (29)  assume 


Ti  =  2^  K  eAx  F(y) 

6j  =  exx  H(y) 

Equations  (44)  and  (45)  then  reduce  to  the  ordinary  differential  equation 
system 


H*  -  Ap2  F  =  4  [XF  +  4  »*  3 

Y-l  a2  M2 

OO  00 

(47) 

F'  +  AH  =  ~  ~  AH 

^  q£ 

with  boundary  conditions 

H ( 0 )  =  H( 1 )  =  0  (48) 
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Eliminating  F  from  equations  (47)  gives  the  single  equation 


H"  +  A2  f*2  H  =  -4t  ^  C0(H"  +  A2  p2  H  +  Xm£  F^  +  o '  ( H '  +  Xm£  F)]  (49) 

oo 

This  equation  is  in  a  form  suitable  for  solution  by  iteration.  The  lowest 
order  solution  (i.e.,  first  approximation)  satisfies 


H 


+  X‘ 


p2  H  =  0 


F 


(50) 


Using  these  results  in  the  RHS  (denoted  by  $)  of  equation  (49)  produces 

H"  +  X2  p2  H  =  -Ar  H'  +  m£  o  H"]  =  (51) 

oo 

The  solution  of  this  equation  provides  an  improved  second  approximation  for 
the  solution  of  equation  (49).  Repetition  of  this  iterative  process  generates 
ever-improving  approximations.  Only  the  solution  for  the  second  approximation 
will  be  given  here.  Whether  or  not  it  is  sufficiently  accurate  can  be 
assessed  by  evaluating  the  numerical  results  presented  later. 

The  solution  for  the  first-approximation  equations  (50)  was  obtained  in 
the  above  isentropic  analysis  as 

oo 

H  =  Z  An  sin  miy  ;  A  =  +  ^  (52) 

1  p 

where  the  An  are  arbitrary  coefficients.  Note  that  the  iterative  solution 
procedure  used  here  is  consistent  in  that  the  isentropic  solution  is  recovered 
for  vanishingly  small  o.  Using  the  approximation  (52)  the  RHS  of  equation 
(51)  can  be  written 

4  00  2  °° 

$  =  — r  [o'  l  n  An  cos  nny  -  n  Mi  0  l  n2  An  sin  nny]  (53) 

Y_i  1  1 
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Equation  (51)  can  be  separated  into  component  equations 


"  2  ? 

H  +  X  p  H  =  4> 
n  n  H  n  Yn 


(n  =  1,  2,...) 


(54) 


where  4>n  represents  the  nth  component  of  the  source  term  $  defined  as 


4  "nAn  r  • 

4>n  =  2-2  ‘-0  cos  ni,y  "  11  o  n  sin  nny] 

^  q» 


(55) 


For  consistency,  the  separation  constant  X must  also  be  expanded  as 


x  =  +  (HI  +  sx  +  . . . ) 
n  -  '  p  n  ’ 


(56) 


so  that  equation  (54)  becomes 


ii  2  ? 

H  +  n  n  H  =  $  -  2npn  6X  A  sin  nny 

n  n  Tn  n  n  J 


(57) 


The  solution  of  this  equation  which  satisfies  the  boundary  conditions  (48)  is 


y 

Hn  =  Kn  sin  nny  +  p  6Xn  An  y  cos  nTtY  +  7^  $  <*>n  sin  ni,(y_ri)  dn 

*  ri 


(58) 


where  Kn  is  an  arbitrary  coefficient  and 


6Xr 


1  f1 

J  ♦„  sin  nnn 


dn 


(59) 


This  solution  is  valid  regardless  of  the  choice  of  sign  in  equation  (56). 

The  solution  for  the  nth  component  of  F  can  be  obtained  from  the  system 
(47)  as 


-Kn  4  0 


Fn  =  t  (1  ~  rrr  -p-p)  cos  ni,Y  “  6An  An  y  sin  nny 


1  ry 

+  — -  \  $n  cos  nn(y-n)  dq] 
nnp  i  - 


(60) 
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The  choice  of  sign  corresponds  to  that  of  equation  (56). 

Using  the  above  results,  the  approximate  solution  of  the  system  (44)  and 
(45)  can  be  written  as 


l  Kin 
1 


-p"  (1  -  -rr  -5%)  cos  nuy 

'  *  1  *£  ,(1B  *  »„>* 

sin  nuy 


+  Z  *2n 
1 


■ZqooMoo  4  a 

— (l  '  —1  X2>  C0S 

P 

oo 

sin  nuy 


e-(°£  +  6Xn)K 


^  f"  o'2?  *  «n).  ♦  ;  ^  f"  e-'2?  *  “n)» 


where 


1  t y 

fn  =  3  6\n  An  y  sin  nuy  -  —  J  <*n  cos  nu(y-n)  dn  (* 

1  y 

hn  =  p  6An  An  y  cos  nuy  +  —  V  <*n  sin  nu(y-q)  dn  U 

3  o 

This  solution  reduces  to  the  isentropic  result  (35)  when  a  is  zero.  Using 
this  solution,  equation  (43)  can  be  integrated  to  obtain  the  entropy 
perturbation  Sj. 

Cascade  Flow 


For  isentropic  cascade  flow,  equations  ( 17 ) - ( 19 )  reduce  to 


_iL_  rill  +  T  Hi]  .  *21  +  T  =  0 

2q00M00  Lax  ay  J  ay  ax 
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(65) 


ael  ,  391  ,  1  rlll  _  ^[l,  _  n 

3x  T  3y  2q  M  ^3y  T  3x  -* 

CO  oo 

If  the  (x,y)  coordinates  are  chosen  such  that  the  flow  is  periodic  in  the  y 
direction  (see  Figure  2),  then  a  separation  of  variables  defined  as 

Ti  =  2<>A  F^>  eini,y 

oo  oo 

e1  =  HU)  e1ni,y  (66) 

£;  =  irnix 

can  be  used.  The  blade  spacing  has  been  taken  as  unity.  Equations  (64)  and 
(65)  are  thereby  transformed  to  the  ordinary  differential  equation  system 

e2( F ’  +  tF)  +  tH'  -  H  =  0 

(67) 

H'  +  tH  -  tF'  +  F  =  0 

Primes  denote  differentiation  with  respect  to  t,. 

The  solution  of  the  system  (67)  has  the  form 

F  =  aeXE>  H  =  eA£;  (68) 

where  the  eigenvalues  are 


(69) 


(70) 


Using  these  results,  the  general  solution  of  the  system  (64)  and  (65)  can  be 
written  as 


16 


Tl 

91 


where 


I  Kln 

oo 

ao 

E  K2n 


P 

1  [ 

2qooM®i  " 
3 
1 


intr(  A^x  +  y) 
e 


inir{A2x  +  y) 


+ 


*1  =  ~ 

T 

a2  i  - 

Note  that  A 


M2  +  iP(i  +  t2) 

oo 

^TT2 

M2  -  ip(l  +  t2) 

00 

32TT2 

and  a2  are  complex  conjugates. 


(71) 


(72) 
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IV.  BOUNDARY  CONDITION  DEVELOPMENT 


Examination  of  equations  ( 1 ) - ( 4 )  shows  that  at  a  subsonic  far  field 
computational  boundary  there  are  three  downstream-running  waves  and  one 
upstream  running  wave.  Therefore,  the  information  available  from  the 
numerical  solution  is  not  complete  and  differs  at  upstream  and  downstream 
boundaries.  The  information  lacking  is  provided  by  the  boundary  conditions. 

If  the  flow  is  supersonic,  all  waves  are  downstream-running  and  specification 
of  boundary  conditions  is  straightforward. 

Far  field  computational  boundary  conditions  (subsonic)  are  developed  in 
this  section  based  on  the  linearized  Euler  solutions  obtained  in  the  previous 
section.  These  solutions  are  assumed  valid  in  the  region  beyond  the  computa¬ 
tional  boundaries  where  nonlinear  effects  are  small.  Within  the  computational 
boundaries  the  full  nonlinear  Euler  equations  must  be  solved  numerically.  The 
boundary  conditions  provide  for  a  smooth  coupling  of  the  nonlinear  and  linear 
solutions  so  that  the  true  conditions  at  infinity  can  be  imposed. 

The  three  specific  cases  analyzed  in  the  previous  section  will  be 
addressed  in  this  section. 


Duct  Flow 


Both  upstream  and  downstream  duct  flow  boundary  conditions  will  be 
derived  for  isentropic  conditions.  The  upstream  boundary  analysis  is  valid 
even  if  there  is  nonlinear  entropy  production  downstream  of  the  boundary 
within  the  computational  domain.  For  non-isentropic  conditions  only 
downstream  boundary  conditions  will  be  derived,  since  non-isentropic 
downstream  flow  is  the  more  common  situation. 

For  isentropic,  constant  pressure  flow  in  the  upstream  or  downstream 
regions,  the  far  field  Mach  number  M^  is  determined  implicitly  from  the  mass 
flow  per  unit  area  w  by  the  relation 
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1  1+i 

M  [1  +  ^  M2]  2  =  w 


The  mass  flow  is  usually  known  (or  can  be  calculated)  for  a  given  duct  flow 
problem.  The  associated  speed  of  sour.d  and  velocity  are 


hi 

M  1+T 

a„  =  (-1) 

w 


boo  “  ^00  Me 


Therefore,  the  far  field  quantities  Q  and  R  appearing  in  expansions  (20)  are 


Qoo  “  boo  +  1  3o 


Roc  —  qoo  ”  i  a0 


The  far  field  entropy  has  been  set  to  Equations  (21)  relate  Q  and  R  to 

the  variable  T. 

For  non-isentropic,  constant  pressure  flow  in  the  downstream  region, 
equations  (37),  (39),  and  (40)  give  the  approximate  far  field  relationships 


Z  -  2  a°° 

Qoo  ~  Poo  1  a  CO  _  1  -  o 


2  -  2  a“ 
R°° =  q°° "  7-T  a°°  ’  y-i  0 


Too  -  3o 


Quadratic  and  higher  terms  in  o  have  been  neglected. 


Isentropic  Conditions  -  At  a  computational  boundary  (assumed  located  at  x=0) 
the  perturbation  flow  variables  can  be  represented  by  the  Fourier  series 
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(77) 


@1  =  l  An  sin  niry 
1 


Ql  =  z  Bn  cos  nny  (78) 

1 


Rl  =  l  Cn  cos  mry  (79) 

1 


Tj  =  I  En  cos  mry  =  l  (Bn-Cn)  cos  n-ny  (80) 

11 

The  boundary  conditions  will  be  developed  from  relationships  between  the 
Fourier  coefficients  An>  Bn,  and  C  .  The  absence  of  modes  corresponding  to 
n=0  for  ,  R^  and  T^  is  related  to  the  fact  that  these  first-order 
perturbations  must  vanish  at  infinity.  Further  discussion  of  this  topic  is 
presented  below  in  conjunction  with  the  boundary  condition  development  for 
cascade  flows. 

The  general  solution  for  linearized  isentropic  flow  is  given  by  equation 
(35).  Applying  this  solution  at  x=0  and  using  the  series  expansions  ( 77 ) - ( 80 ) 
it  follows  that 


A  =  K.  +  K, 
n  In  2n 


2qJ'/L 

B 


(Kin  "  *^2n) 


Therefore, 

*<ln  =  2  +  2q^M00  (®n  "  ^n)] 

*2n  =  2  '  2qSt>D  ^n  ~ 


(81) 


(82) 

(83) 


For  the  region  upstream  of  the  computational  boundary  (i.e.,  x<0),  the 
exponentially  growing  component  of  the  solution  (35)  can  be  suppressed  by 
forcing  K^n  to  be  zero.  This  requires  the  Fourier  coefficients  to  be  related 
by 
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A  =  — J &- 

n  2qooM, 


(Bn  “  Cn) 


(84) 


fhe  linearized  upstream  solution  is  then 


2qooM» 

'Tf 

oo 

=  Z  An 

„  cos  nuy 

p 

_®1_ 

1 

sin  nuy 

emrx/p 


(85) 


For  isentropic  flow  there  are  two  downstream-running  waves  propagating 
information  to  the  upstream  boundary  from  outside  the  computational  domain  and 
one  upstream-running  wave  propagating  information  from  the  numerical  solution. 
Equation  (84)  provides  one  of  the  lacking  pieces  of  information  from  outside 
the  computational  domain;  the  remaining  information  is  provided  by  combining 
equations  (22),  (78)  and  (79)  to  give 

1  ‘  Moo 

Bn  =  i  +  cn  (86) 

Using  the  expansions  (20)  and  the  Fourier  representation  (79),  the 
coefficients  C  are  determined  from 


Rl  =  I  Cn  cos  nuy  =  Rnum  -  Re  (87) 

1 


where  R  is  the  boundary  distribution  of  R  obtained  from  the  nonlinear 
num 

numerical  solution  and  R^  is  given  by  equation  (75).  The  coefficients  An  and 
Bn  are  then  obtained  from  equations  (84)  and  (86).  Using  the  Fourier 
representations  (77)  and  (78),  the  distributions  of  9  and  Q  on  the  boundary 
(i.e.,  the  boundary  conditions)  are  calculated  according  to 


0(3  =  l  An  sin  nuy 
1 


Qb  =  Q„  +  l  Bn  cos  nuy 
1 


(88) 
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For  the  region  downstream  of  the  computational  boundary  (i.e.,  x>0),  the 
exponentially  growing  component  of  the  solution  (35)  can  be  suppressed  by 
forcing  Kj  defined  by  equation  (82)  to  be  zero.  This  requires  the  Fourier 
coefficients  to  be  related  by 


Cn  =  Br 


+  2 3°°^  Ar 


The  linearized  downstream  solution  is  then 


'Tf 

oc 

=  E  An 

-2q„Moc 

_  cos  nny 

p 

01 

1 

sin  nny 

e-nifx/p 


(89) 


(90) 


For  isentropic  flow  there  are  two  downstream-running  waves  propagating 
information  to  the  downstream  boundary  from  the  numerical  solution  and  one 
upstream-running  wave  propagating  information  from  outside  the  computational 
domain.  Equation  (89)  provides  the  lacking  information.  Using  the  expansions 
(20)  and  the  Fourier  representations  (77)  and  (78),  the  coefficients  and 
are  determined  from 


01  E  Ap  sin  n-ny  -  B^um 
1 


Ql  =  l  Bn  cos  nny  = 
1 


Qnum  "  Q® 


(91) 


where  0  and  Q  are  the  boundary  distributions  of  0  and  Q  obtained  from 
num  ^num 

the  nonlinear  numerical  solution  and  Q  is  given  by  equations  (75).  The 

00 

coefficients  Cn  are  then  obtained  from  equation  (89).  Using  the  Fourier 
representation  (79),  the  distribution  of  R  on  the  boundary  (i.e.,  the  boundary 
condition)  is  calculated  according  to 


Rb  =  R®  +  1  cn  cos  n"y 
1 


=  Qnum 


A 

Y-l 


+ 


P  " 


l  An  cos  n-ny 
1 


(92) 
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The  second  form  of  this  relation  requires  only  the  calculation  of  the 
coefficients  . 

Many  numerical  solution  algorithms  for  the  Euler  equations  use  so-called 
characteristic  far  field  boundary  conditions  in  which  6  and  Q  are  specified 

GO  OO 

at  inflow  boundaries  and  R  is  imposed  at  outflow  boundaries.  The  boundary 

oo 

conditions  (88)  and  (92)  therefore  represent  a  logical  first-order  extension 
of  the  widely-used  characteristic  (or  zero-order)  boundary  conditions. 

Non-isentropic  Conditions  -  For  the  case  of  non-isentropic  flow  crossing  a 
downstream  computational  boundary,  the  Fourier  series  representation  (77)  for 
0^  at  the  boundary  is  still  valid  because  6  is  zero  at  the  duct  walls.  The 
variable  o  which  characterizes  the  entropy  distribution  at  infinity  can  be 
represented  by 


o  =  ow  +  I  Dfc  sin(k-l/2)ny  (93) 

1 

where  a  is  the  wall  value  at  y  =  0.  This  Fourier  series  representation 
w 

assumes  an  even  extension  of  a  for  1  <  y  <  2.  The  source  term  component  <$n 
defined  by  equation  (55)  can  then  be  written  as 


2  ^ 

<f>n  =  — —r  5  {  I  (k-1/2)  D|<  [cos(n+k-l/2)ny  +  cos(n-k+l/2)iry] 

P*  <£  k=l 

(94) 

+  n  m£  i  D|<  [cos(n+k-l/2)ny  -  cos(n-k+l/2)-ny]  -  2ow  n  sin  nny} 
k=l 


The  procedure  for  relating  a  to  the  entropy  distribution  at  the  computational 
boundary  and  determining  the  coefficients  will  be  presented  later. 


The  approximate  solution  for  linearized  non-isentropic  flow  is  given  by 

equation  (61).  Since  it  is  only  an  approximate  solution  of  the  approximate 

system  (44)  and  (45),  the  arguments  used  in  the  isentropic  analysis  for 

choosing  K,  and  Kn  to  eliminate  exponentially  growing  solution  components 
In  2n 

cannot  be  used.  That  is,  the  RHS  of  equation  (51)  has  been  treated  as  a  known 
source  term  by  using  a  lower-order  solution  approximation.  As  a  result,  the 
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behavior  of  the  source  term  for  large  x  directly  influences  the  character  of 
the  solution  obtained.  However,  closer  examination  of  the  isentropic  analysis 
shows  that  the  same  results  would  have  been  obtained  had  the  exponentially 
growing  components  simply  been  ignored.  The  remaining  constant  (i.e.,  KJn  or 
K.,n)  associated  with  the  decaying  solution  component  could  have  been 
determined  by  matching  the  solution  at  x=0  with  the  Fourier  series 
representations  on  the  boundary. 


Using  this  approach  and  selecting  only  the  decaying  component  of  the 
solution  (61),  the  solution  for  the  downstream  region  (i.e.,  x>0)  is 


Tl 

CO 

el 

=  |  K„ 

2CJooMoO  u  \ 

— a —  d  -  rrr  7T~2)  cos  n"y 

p  t  1  P^q,*/ 

sin  nny 


e-(^  +  6Xn)x 


+  i 
1 


2qoAo 

B 

hn 


e-(^+«Xn)x 


where 


(95) 


i  y 

fn  =  B  6An  An  y  sin  niry  -  —  J  4>n  cos  nu(y-n)  dn 


(96) 


i  t  y 

hn  =  3  6An  An  y  cos  niry  +  —  j  4»n  sin  ni,(y_n)  dn 

J  0 

6X"  =  J0  ♦"  sin  n,Tri 


(97) 

(98) 


These  last  three  expressions  are  obtained  from  equations  (59),  (62),  and  (63). 
If  the  upstream  region  (i.e.,  x<0)  were  to  be  analyzed  for  non-isentropic 
conditions,  then  the  other  solution  component  would  need  to  be  selected. 
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The  solution  (95)  can  be  matched  to  conditions  at  the  computational 
boundary  x=0  by  letting 


=  An  +  6A 
n  n 


and  using  the  Fourier  series  representation  (77)  to  give 


(99) 


oo  2  1  y 

6A,,,  =  £  —  f  sin  muy  f  $n  sin  nir(n-y)  dn  dy  - 

n=l  3  o  -^o 


°°  1 

2p  £  An  6Xn  1  sin  muy  cos  nny  y  dy 
n=l  ^  o 


(100) 


Evaluation  of  the  integrals  in  this  expression  along  with  those  appearing  in 
equations  (96)-(98)  is  given  in  Appendix  A. 


All  parameters  appearing  in  the  solution  (95)  can  be  evaluated  in  terms 

of  the  Fourier  coefficients  A  and  0.  .  The  A  coefficients  can  be  determined 

n  K  n 

in  the  same  manner  as  for  isentropic  flow  using  equation  (91),  namely. 


oo 

£  An  sin  nny  =  0num  (101) 

1 


where  0  is  the  distribution  of  6  along  the  boundary  obtained  from  the 
nonlinear  numerical  solution.  The  coefficients  can  be  determined  by 
solving  the  entropy  equation  (43).  Neglecting  higher-order  terms,  its  solu¬ 
tion  is 


Sl  *  -  -T-  £  7T  t  !r  sin  n.y  *  *An>x 
y-1  n  dy  .  n 


where,  from  equation  (93), 


jj2  =  IT  £  (k-1/2)  Dk  cos(k-l/2)ny 

ay  i 


(102) 


(103) 
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Combining  these  two  expressions  gives 


Sj  =  -  Z  i  An  Dk  [sin(n+k-l/2)ny  + 

Y  k=l  n=l  n 


sin(n-k+l/2)iry] 


■(r+  6*n)* 


(104) 


Since  entropy  is  convected  downstream  along  streamlines,  the  distribution 
at  the  computational  boundary,  denoted  by  Sn(jm,  is  known  from  the  numerical 
solution.  From  the  expansions  (14)  and  the  definition  (37)  it  follows  that 


S  =  S  +  S, (0 ,y ) 
num  <®  1'  J ' 

(105) 

=  (1  -  o)  +  S^O.y) 

Introducing  the  relationship  (93)  for  o  and  S^O.y)  from  the  solution  (104) 
gives 


I  Dk  sin(k-l/2)ry  =  1  -  ow  -  Snum  - 

k=l  c 

(106) 

^  I  Z  A  D.  [sin(n+k-l/2)uy  +  sin(n-k+l/2)T?y] 

2  k=l  n=l  " 

The  coefficients  must  therefore  satisfy 


Om  =  2  ^  [1  ■  °w  “  Snum]  sin(m*l/2)ny  dy  - 


I  ^  Vk  -  I  =  Dk  #-H-i  *  |  i 

k=l  k=l  k=m+l 


(107) 


This  equation  can  be  solved  for  each  Dm  by  iteration,  started  by  using  the 
first  integral  term  as  an  initial  guess.  Note  that  the  wall  value  is  known 
from  the  numerical  solution  since  the  wall  is  a  streamline. 
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Combining  the  definition  (7)  and  the  expansions  (14)  gives  the 
first-order  relation 

R  =  Q  -  T  -  T.  (108) 

Applying  this  approximation  along  with  the  solution  (95)  at  the  computational 
boundary  x=0,  the  distribution  of  R  on  the  boundary  (i.e.,  the  boundary 
condition)  can  be  calculated  according  to 


4  - 

Bb  =  Onum  “  ^TJ  a°°  + 


2q00M«0 

3 


(1 


-J - 2_) 

Vi  p2qIJ 


(An  +  6An)  cos  nny 


2q<x,M« 

3 


oo 


I  fn 
1 


(109) 


The  approximation  (76)  for  has  been  used  and  Qnum  denotes  the  boundary 
distribution  of  Q  obtained  from  the  nonlinear  numerical  solution.  This 
expression  reduces  to  the  boundary  condition  (92)  if  the  flow  becomes 
isentropic. 


Cascade  Flow 


Both  upstream  and  downstream  cascade  boundary  conditions  will  be  derived 
for  isentropic  conditions.  The  analysis  parallels  that  for  isentropic  duct 
flow  except  that  the  parameters  appearing  in  the  general  linearized  solution 
(71)  are  complex  quantities.  For  isentropic  conditions,  the  cascade  turns  the 
flow  without  losses.  Therefore,  all  upstream  and  downstream  far  field 
quantities  (e.g.,  Q  ,  R^,  etc.)  can  be  determined  from  the  downstream  pressure 
at  infinity,  p^. 


27 


At  a  computation  boundary  (assumed  located  at  x=0)  the  perturbation  flow 
variables  can  be  represented  by  the  general  Fourier  series 


01 


£  An  e 


inuy 


y-  +  £  An  cos  nny  +  £  An  sin  nny 
c  1  1 


Qi  = 


£  Bn  e 


mny 


g  oo  _  oo 

y-  +  £  Bn  cos  nny  +  £  Bn  sin  mry 
£  1  1 


Rl  =  £  Cn  e 


iniry 


y-  +  £  Cn  cos  nny  +  £  Cn  sin  nny 
C  1  1 


(110) 


(111) 


(112) 


Tl  =  £  En  e 


inny 


.  Eo 


=  y  +  £  En  cos  nny  +  £  En  sin  nny 
1  1 


(113) 


The  boundary  conditions  will  be  developed  from  relationships  between  the 

Fourier  coefficients  A  ,  B  and  C  .  F< 

n  n  n 

relationships  (for  n  =  0,  1,  2,  .  .  .) 


Fourier  coefficients  A  ,  B  and  C  .  Fourier  analysis  provides  the  additional 

n  n  n 


Zn  =  Zn  +  Z-n 


Zn  =  ^Zn  '  Z-n> 


Zn  =  \  <Zn  -  1  Zn> 


z-  =  ?  (Zn  +  ’  Z-> 


(114) 
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where  Z  represents  any  of  the  coefficients  A,  B,  C  or  E.  Since  perturbations 

must  vanish  at  infinity,  A  ,  §  ,  C  and  E  must  each  be  zero. 

0  0  0  0 

Applying  the  solution  (71)  at  x=0  and  using  the  series  expansions  (110) 
and  (113)  it  follows  that 


An  =  Kln  +  K2n 


2qooM«.  . 

En  =  x  i  (^2n  "  Kln) 


(115) 
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Therefore , 

Kln  =  \  [An  +  4jC  En] 


K2n  =  2  [A"  "  En-* 


(116) 

(117) 


In  the  region  upstream  of  the  computational  boundary  (i.e.,  x<0)  the 
inflow  angle  e  is  known  so  that  t  defined  by  equation  (16)  is  known.  To 

oo 

suppress  the  exponentially  growing  component  of  the  solution  (71),  the 
coefficients  Kln  and  K2n  must  satisfy 


Kln  =0  (n  >  0) 


K2n  =0  (n  <  0) 


(118) 


This  requires  the  Fourier  coefficients  to  be  related  by 


a  =  -  -ifi- 
"  2q«tt, 


En  =  “  45c  (Bn  '  C,l) 


An  = 


En  =  45c (Bn  "  Cn) 


2qoAo 

The  linearized  upstream  solution  is  then 


(n  >  0) 
(n  <  0) 


(119) 
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The  parameters  X ^  and  X are  defined  by  equations  (72).  The  component 
corresponding  to  n=0  has  been  neglected.  This  will  be  discussed  in  more 
detail  below. 

Equation  (86)  used  for  the  duct  flow  analysis  is  again  applicable. 

Equation  (87)  can  be  generalized  to  determine  the  coefficients  Cn>  That  is, 

C  and  C  are  determined  from 
n  n 

oo  _  oo  _ 

I  Cn  cos  nny  +  Z  Cn  sin  mry  =  Rnum  -  R*  (121) 

1  1 

where  R  is  the  boundary  distribution  of  R  obtained  from  the  nonlinear 
num  J 

numerical  solution.  The  coefficient  C  has  been  neglected.  The  coefficients 

o 

An  and  Bn  are  then  obtained  from  equations  (86)  and  (119)  and  the  general 
Fourier  relationships  (114).  Using  the  Fourier  representations  (110)  and 
(111),  the  distributions  of  0  and  Q  on  the  boundary  (i.e.,  the  boundary 
conditions)  are  calculated  according  to 


@b  =  9®  +  T  An  cos  nny  +  Z  sin  mry 
1  1 


Qb  =  Q«>  +  T  Bn  cos  nny  +  Z  Bn  sin  nny 
1  1 


(122) 


Both  Aq  and  Bq  have  been  set  to  zero  in  the  boundary  conditions  (122).  A  zero 
value  of  Bq  imposed  on  the  numerical  solution  through  the  boundary  condition 
should  tend  to  force  a  zero  value  for  CQ  (and  hence  EQ)  in  the  numerical 
solution  according  to  equation  (86).  It  should  be  noted  from  these 
first-order  relations  that,  within  the  context  of  a  Fourier  representation, 
and  6  can  be  viewed  as  the  mean  values  of  Q  and  0  on  the  computational 

oo 

boundary.  Deviations  from  this  condition  can  be  shown  to  be  a  second-order 
effect.  These  arguments  apply  also  to  duct  flows. 
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For  the  region  downstream  of  the  computational  boundary  (i.e.,  x>0),  the 
outflow  angle  8m  is  unknown  so  that  t  is  also  unknown.  To  suppress  the 
exponentially  growing  component  of  the  solution  (71),  the  coefficients  Kjn  and 
l<2n  must  satisfy 

«ln  =  0  (n  <  0)  K2n  =  °  (n  >  0)  (123) 

This  requires  the  Fourier  coefficients  to  be  related  by 


An  = 


1L 


En  =  - 


_ip 


i£L  p.  =  _l£ 


( Bn  -  Cn)  (n  <  0) 


(124) 


An  = 


2^M«,  En  =  2^MZ  (Bn  '  Cn)  (n  >  0) 
The  linearized  downstream  solution  is  then 
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As  for  the  duct  flow  analysis,  the  coefficients  A  and  B  are  determined 

n  n 

from  the  general  Fourier  representations  (110)  and  (111)  according  to 


oo  oo  — 

I  An  cos  nity  +  E  An  sin  nny  =  0num  -  e*, 

1  1 

(126) 

oo  _  00  _ 

E  Bn  cos  nny  +  E  8n  sin  niry  =  Qnum  -  Q*, 

1  1 

where  e  and  Q  are  the  boundary  distributions  of  0  and  Q  obtained  from 
num  xnum  J  _ 

the  nonlinear  numerical  solution.  The  coefficients  A  and  B  have  been 

o  o 

neglected.  The  outflow  angle  0^  can  be  calculated  from 
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(127) 


e»  =  (  ®num  dy 
J  o 

which  is  the  mean  value  of  9  on  the  boundary.  This  justifies  setting  the 

coefficient  A  to  zero.  The  coefficients  C  are  obtained  from  equations 
o  n 

(124).  Using  the  Fourier  representation  (112),  the  distribution  of  R  on  the 
boundary  (i.e.,  the  boundary  condition)  is  calculated  from 


Rb  =  Roo  +  E  Cn  cos  nny  +  E  Cn  sin  n-ny 
1  1 

(128) 

4  2qMM<»  rw  ~  00  -  . 

=  Qnum  ■  rrr  a-  +  — 5 —  [E  An  cos  n-ny  -  E  An  sin  nny] 

Y  1  P  1  1 

The  second  form  of  this  relation  requires  only  the  evalution  of  the 
coefficients  An  and  A^.  It  is  analogous  to  equation  (92)  for  duct  flow. 

The  coefficient  CQ  has  been  set  to  zero  which,  when  imposed  on  the 

numerical  solution  through  the  boundary  condition,  should  tend  to  force  a  zero 

value  for  B  (and  hence  E  ).  That  is,  the  boundary  condition  tends  to  enforce 
o  0 

R  as  the  mean  value  of  R  on  the  boundary.  This  argument  applies  also  to  duct 

00 

flows. 
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V.  APPLICATIONS 


Numerical  solutions  of  the  Euler  equations  have  been  calculated  for 
two-dimensional,  steady  duct  flows  using  the  first-order  boundary  condition 
procedures  developed  in  the  previous  section.  Both  isentropic  and 
non-isentropic  cases  were  analyzed.  Non-isentropic  flow  was  produced  by  a 
shock  wave  located  in  a  nozzle  portion  of  the  duct.  First-order  boundary 
condition  results  are  compared  with  those  produced  using  the  conventional 
zero-order  characteristic  boundary  conditions. 

A  solution  algorithm  was  used  for  the  nonlinear  Euler  equations  (l)-(4) 
which  is  based  on  the  method  presented  in  Reference  1.  It  uses  explicit  time 
integration  to  relax  to  steady  state  conditions  and  includes  the  shock-fitting 
procedure  of  Reference  4  to  accurately  calculate  flows  containing  shock  waves. 
She  k  waves  appear  in  the  solution  as  distinct  discontinuities  which  satisfy 
the  Rankine-Hugoniot  relations.  It  should  be  noted  that  the  boundary 
condition  analysis  is  independent  of  the  choice  of  inviscid,  nonlinear 
solution  method. 

The  duct/nozzle  geometry  is  shown  schematically  in  Figure  1.  The  flow  is 
characterized  by  p  ,  the  downstream  pressure  at  infinity,  which  produces  a 

oo 

mass  flow  per  unit  area  w  through  the  duct.  The  linearized  solutions  given  by 
equations  (85),  (90)  and  (95)  are  assumed  valid  in  the  semi-infinite  regions  I 
and  III  and  the  computational  boundary  conditions  are  applied  at  the  upstream 
and  downstream  boundaries  AA  and  BB  of  the  nonlinear  computational  region  II. 

The  actual  shape  of  the  duct/nozzle  and  the  computational  grid  are  shown 
in  Figure  3.  The  nozzle  contour  is  sinusoidal  and  symmetric  about  the 
centerline.  The  computational  grid  for  this  portion  of  the  nozzle  had 
dimensions  41  x  21.  The  area  ratio  is  .75  and  the  upstream  and  downstream 
areas  are  equal.  For  these  constant  area  sections  of  the  duct,  additional 
rectangular  grid  cells  could  be  added  without  altering  the  basic  41  x  21  grid. 
This  served  to  minimize  the  effect  of  grid  changes  on  the  calculations  when 
the  computational  boundaries  were  moved  in  order  to  assess  the  accuracy  of  the 
boundary  conditions. 
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8ecause  the  configuration  is  symmetric,  calculations  were  limited  to  the 
lower  half  of  the  duct  and  a  centerline  symmetry  condition  was  used.  Although 
the  configuration  used  for  these  calculations  is  simple,  the  boundary 
condition  analysis  of  the  previous  section  is  general  and  applicable  to 
unsymmetric  configurations  having  unequal  upstream  and  downstream  areas.  Use 
of  the  simple  configuration  is  sufficient  to  demonstrate  the  validity  of  the 
analysis. 


Isentropic  Conditions 

For  isentropic  flow  the  first-order  upstream  and  downstream  boundary 
conditions  are  given  by  equations  (88)  and  (92),  respectively.  The  associated 
analytic  far  field  solutions  are  given  by  equations  (85)  and  (90).  The 
zero-order  (or  characteristic)  boundary  conditions  consist  of  imposing  the 
constant  value  of  Qro  and  a  zero  value  of  9  along  the  upstream  boundary  and  the 
constant  value  of  R  along  the  downstream  boundary.  The  values  of  Q  and  R 
are  determined  from  equations  (74)  and  (75). 

Results  are  presented  for  a  single  value  of  p^  but  with  the  computational 
boundaries  located  at  several  different  longitudinal  stations.  The  relative 
accuracy  of  the  zero  and  first-order  boundary  conditions  can  then  be 
evaluated. 

Case  1  -  Results  obtained  using  the  complete  grid  shown  in  Figure  3  are 
presented  in  Figures  4  and  5.  This  grid  has  20  columns  of  grid  cells  in  both 
the  upstream  and  downstream  constant  area  portions  of  the  duct.  For  this  case 
the  computational  boundaries  were  far  enough  upstream  and  downstream  that  the 
zero-order  and  first-order  boundary  condition  results  were  nearly  identical. 
Figure  4  shows  pressure  and  Mach  number  distributions  along  the  centerline  and 
lower  wall  of  the  duct/nozzle.  Pressure,  Mach  number,  and  flow  angle  contours 
are  presented  in  Figure  5.  These  results  serve  as  a  reference  for  evaluating 
the  accuracy  of  solutions  where  the  computational  boundaries  are  moved  closer 
to  the  nozzle  portion  of  the  duct. 

Case  2  -  Results  for  a  shortened  duct  obtained  using  the  zero-order  boundary 
conditions  are  presented  in  Figure  6  and  7.  There  were  5  columns  of  grid 
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cells  in  both  the  upstream  and  downstream  constant  area  portions  of  the  duct 
for  this  case.  Pressure  and  Mach  number  distributions  are  shown  in  Figure  6 
and  contours  are  shown  in  Figure  7.  The  most  noticeable  effect  of  the 
boundary  proximity  is  a  change  in  peak  values  of  pressure  and  Mach  number  and 
a  significant  amount  of  longitudinal  asymmetry. 

Case  3  -  The  previous  shortened  duct  case  was  recalculated  using  the 
first-order  boundary  conditions.  These  results  are  presented  in  Figures  8  and 
9.  Pressure  and  Mach  number  distributions  are  shown  in  Figure  8  and  contours 
are  shown  in  Figure  9.  The  results  within  the  numerical  solution  portion  of 
the  domain  are  almost  identical  to  those  shown  in  Figures  4  and  5.  Linearized 
solution  results  obtained  from  equations  (85)  and  (90)  have  been  added 
upstream  and  downstream  of  the  computational  boundaries.  It  is  evident  that 
the  linearized  far  field  analytic  solutions  provide  for  a  smooth  transition 
across  the  computational  boundary  to  the  true  far  field  conditions  at 
infinity. 

Case  4  -  Results  for  a  more  drastic  case  using  the  first-order  boundary 
conditions  are  presented  in  Figures  10  and  11.  Only  one  column  of  grid  cells 
was  used  in  the  upstream  and  downstream  constant-area  portions  of  the  duct. 
Pressure  and  Mach  number  distributions  are  shown  in  Figure  10  and  contours  in 
Figure  11.  Even  with  the  computational  boundaries  extremely  close,  the 
numerical  solution  in  the  nozzle  portion  of  the  duct  is  affected  very  little. 

Boundary  Comparison  -  A  more  quantitative  comparison  between  the  zero-  and 
first-order  boundary  conditions  can  be  obtained  by  examining  the  distribution 
of  flow  variables  along  the  computational  boundaries.  Figure  12  shows  a 
comparison  of  the  distribution  of  pressure  and  flow  angle  along  the  upstream 
and  downstream  boundaries  for  Cases  2  and  3.  Results  taken  from  the  numerical 
solution  of  Case  1  at  the  same  longitudinal  locations  are  also  shown.  There 
is  little  difference  between  the  Case  1  and  Case  3  results.  The  larger 
differences  between  the  Case  2  and  Case  3  pressures  occurring  along  the 
upstream  boundary  are  primarily  due  to  the  fact  that  a  zero  flow  angle  is 
imposed  along  the  boundary. 
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Implementation  of  the  first-order  boundary  conditions  obviously 
introduces  additional  computational  effort.  However,  the  boundary  conditions 
do  not  have  to  be  updated  with  every  iteration  of  the  numerical  solution.  The 
small  additional  effort  is  more  than  offset  by  the  large  reduction  in  number 
of  grid  points  required  and  the  fewer  solution  iterations  necessary  for 
convergence  because  information  propagates  between  the  computational 
boundaries  more  quickly.  For  the  isentropic  results  presented  above,  only  5 
Fourier  modes  were  required  to  give  a  reasonably  accurate  representation  of 
the  boundary  distributions  of  R  (upstream)  and  0  (downstream)  for  the  Case  3 
calculations.  For  the  Case  1  and  Case  4  calculations,  the  number  of  Fourier 
modes  required  were  3  and  7,  respectively. 

Non-Isentropic  Conditions 

Non-isentropic  flow  was  produced  by  lowering  the  value  of  p^  so  that  a 
shock  wave  formed  in  the  nozzle  portion  of  the  duct.  For  this  situation  the 
isentropic,  first-order  upstream  boundary  conditions  (88)  and  the  associated 
linearized  solution  can  be  used.  The  upstream  far  field  quantities  Q  and  R 
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are  constant  and  must  be  determined  from  the  mass  flow  through  the  nozzle  as 
explained  when  deriving  equations  (75).  The  mass  flow  can  be  determined  by 
numerically  integrating  the  mass  flux  crossing  the  downstream  computational 
boundary  each  time  the  boundary  conditions  are  updated. 

The  first-order  downstream  boundary  condition  is  given  by  equation  (109) 
and  the  associated  analytic  solution  by  equation  (95).  It  should  be  noted 
that  this  solution  is  only  an  approximation  to  the  solution  of  the  linearized 
equations  (23)-(25).  The  downstream  far  field  quantities  and  vary  in 
the  streamline-normal  direction  and  can  be  determined  from  equations  (76). 
Zero-order  downstream  boundary  conditions  as  used  herein  consist  of  imposing  a 
distribution  of  R^  along  the  boundary  given  by  equations  (76)  with  o  replaced 
by  1-  S  .  This  in  itself  represents  an  improvement  over  some 
applications  of  zero-order  characteristic  boundary  conditions  in  which  a 
constant  value  of  the  Riemann  variable  associated  with  upstream-running  waves 
is  imposed  along  the  boundary. 

Results  are  presented  for  two  different  values  of  p^  chosen  such  that  a 
weak  and  a  strong  shock  are  produced  in  the  nozzle.  The  strong  shock  spans 
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the  width  of  the  nozzle  while  the  weak  shock  does  not.  Two  computational 
boundary  locations  are  used  for  each  value  of  p^  so  that  the  accuracy  of  the 
non-isentropic  boundary  conditions  can  be  evaluated. 

Case  5  -  Weak  shock  results  obtained  using  the  complete  grid  of  Figure  3  are 
shown  in  Figures  13  and  14.  Pressure  and  Mach  number  distributions  are  shown 
in  Figure  13  for  the  centerline  and  lower  wall.  Pressure,  Mach  number, 
entropy  and  flow  angle  contours  are  presented  in  Figure  14.  The  non-constant 
Mach  number  in  the  downstream  far  field  is  evident  in  Figure  13.  Note  that 
the  scales  used  in  Figure  13  are  much  more  compressed  than  those  used  for  the 
isentropic  results. 

For  this  case,  there  were  slight  differences  produced  by  the  choice  of 
boundary  conditions,  so  these  two  sets  of  results  are  shown  separately.  These 
differences  can  be  attributed  to  difficulties  with  the  numerical  fitting  of  a 
weak  shock  discontinuity.  The  shock  fitting  procedure  is  extremely  sensitive 
to  small  disturbances  in  the  numerical  solution  when  the  shock  is  weak, 
especially  in  the  vicinity  of  the  tail  of  the  shock  located  near  the  center  of 
the  computational  domain.  These  interactions  were  stronger  when  the 
zero-order  boundary  conditions  were  used. 

Placement  of  the  computational  boundary  even  further  downstream  by  using 
additional  columns  of  grid  cells  caused  very  little  change  to  the  solution 
when  the  first-order  boundary  conditions  were  used,  but  produced  noticeable 
differences  when  the  zero-order  conditions  were  used.  For  this  reason  the 
results  of  Figures  13b  and  14b  can  be  considered  more  accurate  and  can  serve 
as  a  reference  for  evaluating  the  accuracy  of  solutions  where  the 
computational  boundaries  are  closer  to  the  nozzle  portion  of  the  duct. 

Case  6  -  Results  for  a  shortened  duct  using  the  zero-order  boundary  conditions 
are  presented  in  Figures  15  and  16.  There  were  5  columns  of  grid  cells  in 
both  the  upstream  and  downstream  constant  area  portions  of  the  duct  for  this 
case.  Pressure  and  Mach  number  distributions  are  shown  in  Figure  15  and 
contours  are  shown  in  Figure  16.  The  effect  of  the  boundary  proximity  is 
clearly  evident  in  these  results. 
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Case  7  -  The  short  duct  flowfield  was  recalculated  using  the  first-order 
boundary  conditions.  These  results  are  presented  in  Figures  17  and  18. 
Pressure  and  Mach  number  distributions  are  shown  in  Figure  17  and  contours  in 
Figure  18.  The  results  agree  closely  with  those  shown  in  Figures  13b  and  14b 
for  the  numerical  solution  portion  of  the  domain.  Linearized  solution  results 
obtained  from  equations  (85)  and  (95)  have  been  added  upstream  and  downstream 
of  the  computational  boundaries.  Even  though  the  solution  (95)  is  appro  '- 
mate,  having  been  obtained  by  linearizing  the  thermodynamic  relation  (36)  and 
using  only  a  two  iteration  approximation  for  the  solution  of  equation  (49),  it 
still  leads  to  a  useful  boundary  condition  and  an  adequate  prediction  of  the 
downstream  flowfield. 

Case  8  -  Strong  shock  results  obtained  using  the  complete  grid  of  Figure  3  are 
presented  in  Figures  19  and  20.  Pressure  and  Mach  number  distributions  are 
shown  in  Figure  19  and  contours  in  Figure  20.  The  significant  variation  of 
Mach  number  in  the  downstream  portion  of  the  duct  is  clearly  evident  in  Figure 
19.  The  shock  fitting  procedure  was  much  less  sensitive  to  numerical 
disturbances  for  this  strong  shock  case.  As  a  result,  the  solutions  obtained 
using  the  zero-order  and  first-order  boundary  conditions  for  the  long  duct 
were  nearly  identical.  These  results  can  serve  as  a  reference  for  evaluating 
the  accuracy  of  solutions  where  the  computational  boundaries  are  closer  to  the 
nozzle  portion  of  the  duct. 

Case  9  -  Results  for  the  shortened  duct  using  the  zero-order  boundary 
conditions  are  presented  in  Figures  21  and  22.  Pressure  and  Mach  number 
distributions  are  shown  in  Figure  21  and  contours  in  Figure  22.  The 
inaccuracy  of  the  zero-order  boundary  conditions  for  this  case  is  evident. 

One  noticeable  effect  is  the  change  in  shock  position. 

Case  10  -  The  short  duct  flowfield  was  recalculated  using  the  first-order 
boundary  conditions.  These  results  are  presented  in  Figures  23  and  24. 
Pressure  and  Mach  number  distributions  are  shown  in  Figure  23  and  contours  in 
Figure  24.  For  the  numerical  solution  portion  of  the  domain,  the  results 
agree  closely  with  those  shown  in  Figures  19  and  20.  Linearized  solution 
results  obtained  from  Equations  (85)  and  (95)  have  been  added  upstream  and 
downstream  of  the  computational  boundaries.  The  transition  across  the 
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boundary  is  reasonably  smooth  indicating  that  the  approximations  used  in 
developing  the  non-isentropic  boundary  conditions  do  not  produce  a  strong 
mismatch  between  the  numerical  and  analytic  solutions. 

Boundary  Comparison  -  A  more  quantitative  comparison  between  the  zero-  and 
first-order  downstream  boundary  conditions  can  be  obtained  from  the 
distribution  of  flow  variables  along,  the  boundary.  Figure  25  shows  a 
comparison  of  the  pressure  and  flow  angle  distributions  along  the  downstream 
boundary  for  Cases  6  and  7  (weak  shock)  and  Cases  9  and  10  (strong  shock). 
Results  taken  from  the  numerical  solutions  of  Case  5  (first-order  boundary 
conditions)  and  Case  8  at  the  same  longitudinal  location  are  also  shown.  The 
validity  of  the  first-order  boundary  condition  analysis  is  clearly  evident. 

The  minor  deviations  are  due  in  part  to  the  approximations  used  in  obtaining 
the  solution  (95). 

For  the  non-isentropic  calculations,  9  Fourier  modes  were  required  to 
give  a  reasonably  accurate  representation  of  the  downstream  boundary 
distributions  of  0  and  a  for  the  Case  7  and  10  calculations.  For  the  Case  5 
and  8  calculations,  only  7  modes  were  used.  Implementation  of  the  first-order 
boundary  conditions  for  non-isentropic  situations  clearly  requires  more 
computational  effort  than  for  isentropic  flow.  However,  this  increase  is 
warranted  in  view  of  the  fact  that  the  far  field  flow  can  contain  significant 
streamline-normal  variations  in  quantities  other  than  pressure  (e.g.,  see 
Figure  19). 
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VI.  SUMMARY 


Far  field  computational  boundary  conditions  have  been  developed  for  2D 
duct  and  cascade  flows.  These  first-order  boundary  conditions  are  derived 
from  analytic  solutions  of  the  linearized  Euler  equations  and  represent  a 
logical  extension  of  the  zero-order  (or  characteristic)  boundary  conditions 
commonly  used  in  the  numerical  solution  of  nonlinear  fluid  dynamic  equations. 
The  boundary  conditions  and  analytic  solutions  provide  a  smooth  transition 
across  a  computational  boundary  to  the  true  far  field  conditions  at  infinity. 
The  boundary  procedure  is  general  in  that  it  can  be  used  in  conjunction  with 
any  numerical  solution  method. 

For  the  case  of  isentropic  flow  the  linearized  Euler  equations  are  solved 
exactly  for  duct  and  cascade  flow  using  separation  of  variables  and  Fourier 
analysis.  For  non-isentropic  flow  the  far  field  can  contain  significant 
streamline-normal  variations  in  quantities  other  than  pressure  and  flow  angle, 
and  these  variations  cannot  be  treated  as  small  perturbations  from  uniform 
flow  conditions.  For  the  case  of  duct  flow  the  non-isentropic  linearized 
equations  are  solved  approximately  using  an  approximate  thermodynamic  relation 
and  an  iterative  technique.  The  validity  of  this  approximate  solution  has 
been  verified  by  numerical  results.  Use  of  non-constant  zero-order  boundary 
conditions,  which  vary  on  the  boundary  as  a  result  of  entropy  variations,  also 
offer  substantial  improvement  over  constant  conditions  derived  solely  from  the 
far  field  pressure.  Extension  of  the  non-isentropic  analysis  to  cascade  flows 
is  straightforward. 

Use  of  zero-order  (or  characteristic)  boundary  conditions  requires  that 
the  computational  boundaries  be  located  far  from  the  nonlinear  region  of  the 
flow.  Closer  placement  of  the  boundaries  may  result  in  a  significant  amount 
of  solution  degradation.  The  first-order  boundary  conditions  allow  the 
boundaries  to  be  located  much  closer  thereby  reducing  the  number  of  grid 
points  needed  for  the  numerical  solution  and  also  the  number  of  iterations  for 
solution  convergence.  This  allows  a  significant  reduction  in  the  amount  of 
computational  effort  required  for  the  nonlinear  numerical  solution  because  the 
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additional  calculations  required  for  the  first-order  boundary  conditions  is 
modest. 


The  boundary  condition  procedures  developed  herein  can  be  extended  to 
axisymmetric,  three-dimensional,  and  viscous  flows,  including  external  flows. 
It  is  recommended  that  these  analyses  be  undertaken  and  also  that  the  present 
results  be  extended  to  non-isentropic  cascade  flows. 
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III.  APPENDIX  A:  EVALUATION  OF  INTEGRALS 


Evaluation  of  the  integrals  appearing  in  equations  (96),  (97),  (98)  and 
(100)  will  be  presented  here.  They  are  defined  as  shown  below. 

Equation  (96): 

i  y 

II  =  rr  C  $n  cos  nu(y-n)  dn 
3  o 

Equation  (97): 

1  y 

*2  =  rvn  (  *n  s1n  dr| 

3  o 

Equation  (98): 

‘3  =  *n  n,y  dy 

Equation  (100): 

1 4  =  f  sin  mny  cos  nny  y  dy 

3  o 

1  y 

15  =  f  sin  miry  f  $n  sin  nn(n-y)  dn  dy 

^  0  '0 

The  quantity  $n  appearing  in  the  integrands  is  defined  by  equation  (94). 
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Evaluation  of  Ij 


*1  =FI 


P2qi 


{  £  (k 
k=l 


1 .  n  rsin(n+k-l/2)Tiy  +  sin  rnt.v 
Uk  L  2n+k-l/2 


,  sin(n-k+l/2)uy  +  sin  niry  ,  sin(n+k-l/2hy  -  sin{n-k+l/2)irv-i 
2n'k+l/2  k- 1/2  J 

,  _  „2  “  rv  rsin(n+k-l/2)iry  +  sin  niry  si n ( n-k+1/2 ) -n v  +  sin  niry 

n  Uk  L  2n+k-l/2  '  2n-k+l/2 

,  sin(n+k-l/2)ny  +  sin(n-k+l/2)iry  -  2sin  n-ny-i  „.-.2  .  •. 

+  - 1 - -L-JL - v7T/2 — L~J~JL - *3  "  2M~  °w  niry  sin  niry} 


Evaluation  of  I2 

12  =  -4  -fe?  (  7  rk  -  4  ot  r--^n— 1/-2)it/--co^  nm 

P2ql  lk=i  2  k  L  2n+k-l/2 

cos(n-k+l/2)Tiy  -  cos  ntiy  cos(n-k+l/2)ny  -  cos(n+k-l/2)ny  , 
2n-k+l/2  k- 1/2  J 

"2  “  n  rcos(n+k-l/2)TT.y  -  cos  nny  cos(n-k+l/2)Tiy  -  cos  niry 

n  k=l  L  2n+k-l/2  '  2n-k+l/2 

cos(n+k-l/2)iry  +  cos ( n-k+1/2 )iry  -  2cos  niry 

k- 1/2  J 

-2 

-  2M*.  ow  [sin  nity  -  nny  cos  nny]} 

Evaluation  of  I3 

1  1  °° 

13  "  vT  3^2  {k4(k  ’  2*  °k  4n+k-l/2  +  2n- k+1/2^ 

oo  „ 

+  °  ^  k  =  l  °k  4n+k-l/2  '  '2n-k+l/2  '  kTT72^  n  °w} 
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Evaluation  of  I4 


,  1 4m+n 

(li)  m 

11  m2-n2 


Evaluation  of  I5 


-  |  p  Am  6Am  -  |  J1  (1-y)  <t»m  cos  rrnty  dy  ;  n=m 

15  = 

•y  y  f1  $n  sin  mi,y  dv  ■  (“l)m+n  P  rn  An  6Xn]  ;  n*n 
J0 

1  1  m  Ani  11  1 

(1-y)  cos  my  dy  =  — r  {  1  (k-^)  Du  [ - ~  +  - - — 5  + 

0  m  ^_1  p2q2  k=i  2  *  (2m+k-l/2)2  (2m-k+l/2)2 

OO  r'  x 

- — x]  +  m  M2  I  Dk  [ - - - X - - - 5]  -  n  ow  m£) 

(k-1/2)2  00  k=1  <2m+k-l/2)2  (2m-k+l/2)2 

1  #n  sin  m.y  dy  =  tlrEt<k“5)  °k  [m+nlk-i72  +  m-dU/2  +  m+n-k+T/2 

O  P  Hoc  k=l 

+  m-n+k-l/2]  +  "  ^  ^  Dk  ^m+n+k-1/2  +  m-n-k+1/2 


m+n-k+1/2  m-n+k-1/2- 
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Figure  1.  Duct/Nozzle  Schematic 


Figure  2.  Cascade  Schematic 


Figure  3.  Duct/Nozzie  Geometry  and  Grid 


46 


ci«..ro  d  Pressure  and  Mach  Number  Distributions 
Figure  4.  First. order  Boundary  Conditions 

^  _  #1  AA 


p  =0.90 

'oo 


47 


Computational  Grid  (81x21) 
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Figure  5.  Pressure,  Mach  Number  and  Flow  Angle  Contours 
Zero-  and  First-Order  Boundary  Conditions 
p  =  0.90 
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Computational  Grid  (51x21) 


Figure  7.  Pressure,  Mach  Number  and  Flow  Angle  Contours 

Zero-Order  Boundary  Conditions 
p  =  0.90 
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Mach  Number  Pressure 


Mach  Number  Contours 


Computational  Grid  (43x21) 


Figure  11.  Pressure,  Mach  Number  and  Flow  Angle  Contours 

First-Order  Boundary  Conditions 
p  =  0.90 
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Mach  Number  Pressure 


Figure  13a.  Pressure  and  Mach  Number  Distributions 

Zero-Order  Boundary  Conditions 

p  =0.83 
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Figure  13b.  Pressure  and  Mach  Number  Distributions 
First-Order  Boundary  Conditions 
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Figure  14a.  Pressure,  Mach  Number,  Entropy  and 
Flow  Angle  Contours 
Zero-Order  Boundary  Conditions 
p  =0.83 
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Pressure  Contours 


Mach  Number  Contours 


Entropy  Contours 


Flow  Angle  Contours 


Computational  Grid  (81x21) 
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Figure  14b.  Pressure,  Mach  Number,  Entropy  and 

Flow  Angle  Contours 

First-Order  Boundary  Conditions 

p  =  0.83 
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Pressure  Contours 


Mach  Number  Contours 


Entropy  Contours 


Flow  Angle  Contours 


Computational  Grid  (51x21) 


Figure  16.  Pressure,  Mach  Number,  Entropy  and 
Flow  Angle  Contours 
Zero-Order  Boundary  Conditions 
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Computational  Grid  (51x21) 


Figure  17.  Pressure  and  Mach  Number  Distributions 
First-Order  Boundary  Conditions 
p  =  0.83 
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Pressure  Contours 


Mach  Number  Contours 


Entropy  Contours 


Computational  Grid  (51x21) 


Figure  18.  Pressure,  Mach  Number  ntropy  and 
Flow  Angle  Contours 
First-Order  Boundary  Conditions 
p  =0.83 
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Figure  19.  Pressure  and  Mach  Number  Distributions 
Zero-  and  First-Order  Boundary  Conditions 
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Computational  Grid  (81x21) 


Figure  20.  Pressure,  Mach  Number,  Entropy  and 
Flow  Angle  Contours 

Zero-  and  First-Order  Boundary  Conditions 
p  =  0.78 


Figure  21.  Pressure  and  Mach  Number  Distributions 
Zero-Order  Boundary  Conditions 
p  =0.78 
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Pressure  Contours 


Entropy  Contours 


Flow  Angle  Contours 


Computational  Grid  (51x21) 


Figure  22.  Pressure,  Mach  Number,  Entropy  and 
Flow  Angle  Contours 
Zero-Order  Boundary  Conditions 
p  =0.78 
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Mach  Number  Pressure 


Computational  Grid  (51x21) 


Figure  23.  Pressure  and  Mach  Number  Distributions 
First-Order  Boundary  Conditions 
p  =0.78 
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Pressure  Contours 


Mach  Number  Contours 


Entropy  Contours 


Flow  Angle  Contours 


Computational  Grid  (51x21) 


Figure  24.  Pressure,  Mach  Number,  Entropy  and 
Flow  Angle  Contours 
First-Order  Boundary  Conditions 
p  =  0.78 
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Pressure  Distribution 
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Figure  25.  Pressure  and  Flow  Angle  Distributions 
on  Downstream  Computational  Boundary 
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